Anomalous electronic structure and pseudogap effects in Ndi.gsCeo.isCuC^ 
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We report a high-resolution angle-resolved photoemission spectroscopic (ARPES) study of the 
electron-doped (n-type) cuprate superconductor Ndi.gsCeo.isCuCU. We observe regions along the 
Fermi surface where the near-i?F intensity is suppressed and the spectral features are broad in a man- 
ner reminiscent of the high-energy "pseudogap" in the underdoped p-type (hole doped) cuprates. 
However, instead of occurring near the (it, 0) region, as in the p-type materials, this pseudogap 
falls near the intersection of the underlying Fermi surface with the antiferromagnetic Brillouin zone 
boundary. 

PACS numbers: 79.60.Bm, 73.20.Dx, 74.72.-h 



The electron-doped cuprate superconductors provide a 
unique opportunity to study the physics of doped Mott 
insulators and high-T c superconductors. In addition to 
possessing interesting physics in their own right, the n- 
type materials, with their different normal state proper- 
ties and phase diagram, offer an alternative venue to test 
various theories of high-T c superconductivity. 

Nd2- a; Cei C Cu04± ( 5 is a member of the small family of 
cuprate superconductors that can be doped with elec- 
trons pi. Only an approximate symmetry in the phase 
diagram exists about the zero doping line between p- and 
n-type, as the antiferromagnetic phase is much more ro- 
bust in the electron-doped material and persists to much 
higher doping levels. Superconductivity occurs in a dop- 
ing range that is almost five times narrower. In addition, 
these two ground states occur in much closer proxim- 
ity to each other. Experiments show other contrasting 
behavior between n-type superconductors and their p- 
type counterparts, for instance a T 2 dependence of the 
in-plane resistivity (2J , the lack of incommensurate neu- 
tron scattering peaks H and a temperature induced sign 
reversal of the Hall coefficient [Q. Although they ap- 
pear to share a superconducting d-wave pairing symme- 
try H||], the 2A sc /fcsT c ratio of the n-type materials 
is much smaller than that of its optimally doped p-type 
counterparts J7|. 

Here we report a detailed angle resolved photoemission 
spectroscopic (ARPES) study of the electronic structure 
of this electron-doped cuprate superconductor. Signifi- 
cant differences are found between it and the p-type ma- 
terials that shed light on a number of important topics 
in the high-T c superconductors. We find that the Lut- 
tinger volumed Fermi surface is truncated into several 
pieces with a high-energy "pseudogap" -like suppression 



that forms, not at the maximum of the d-wave functional 
form as in the underdoped p-type materials, but near the 
intersections of the Fermi surface with the antiferromag- 
netic Brillouin zone (AFBZ) boundary. 

Single crystals of Ndi.gsCeo.isCuO^ (NCCO) were 
grown by the traveling solvent floating zone method in 
4 ATM of 02- Details of this growth can be found else- 
where [g. The resulting crystals show an onset of su- 
perconductivity at 25K and a superconducting volume 
(Meissner shielding) of almost 100% at 20K. A second 
batch of crystals, grown at Stanford University under 
similar conditions, showed an onset at 24K with simi- 
larly narrow transition widths. The photoemission data 
obtained between the two batches are identical. 

As it is imperative to determine which characteristics 
of the data are intrinsic and which are possible artifacts 
due to, for instance, the matrix elements for photoexcita- 
tion, data were collected on two different photoemission 
systems. On beamline 5-4 at the Stanford Synchrotron 
Radiation Lab (SSRL), data were taken at 16.5 eV pho- 
ton energy (Tuj) with an incident angle of approximately 
45 degrees and the in-plane polarization along the Cu-0 
bonds (Energy resolution AE ~10 meV and angular res- 
olution A0 ~ 0.5°). At beamline 10.0.1 at the Advanced 
Light Source (ALS) an incident energy of 55eV was used 
in a glancing incidence geometry with in-plane polariza- 
tion at 45 degrees to the Cu-0 bonds (AE ~ 20meV, A0 
~ 0.25°). All displayed spectra were taken at low tem- 
perature (10 — 20K ) in the superconducting state; this 
does not change the conclusions we make regarding the 
normal state electronic structure as the changes in spec- 
tra with the onset of superconductivity are very subtle 
in NCCO [pi. The chamber pressure was lower than 4 x 
10~ n torr. Cleaving the samples at low temperature in 



situ results in shiny flat surfaces which low-energy elec- 
tron diffraction (LEED) [Fig. 2(c)] reveal clean and well 
ordered with a symmetry commensurate with the bulk. 
No signs of surface aging were seen for the duration of 
the experiment (~ 24 hours). 

In Fig.'s 1(a) - (c) we display energy distribution curves 
(EDCs) from high symmetry lines in the Brillouin zone 
(BZ) taken on the SSRL apparatus. In Fig. 1(a), one sees 
a dispersion along the T to (71", it) direction which is ubiq- 
uitous among the cuprates. A broad feature disperses 
quickly towards the Fermi energy (Ep), sharpens to a 
sharp peak at kp (0.467r,0.467r), and then disappears. 
Along the T to (it, 0) direction, the n-type spectra are 
quite different from their p-type counterparts. While the 
low-energy feature in the optimally doped p-type com- 
pounds disperses quite close to Ep and forms a flat band 
region with a correspondingly high density of states close 
to Ep , in NCCO this flat band region is at about 300 meV 
higher energy ||. Along (71", 0) to (71-, it) [Fig. 1(c)] the 
peak disperses to Ep and gives a small peak at kp. The 
spectra also possess a large background contribution that 
has a weak maximum at 300 meV, as seen in the 100% 
curves in Figs. 1(a) and 1(c). 

In addition to the spectra from the high-symmetry di- 
rections, more comprehensive intensity maps were ob- 
tained over a large region in momentum space at both 
photon energies. In Fig. 2(a) and 2(b) we have plotted 
the integrated spectral weight of the EDCs (Hlu — 16.5 
and 55eV respectively) from within a 30 meV window 
about Ep as a function of k. This gives a measure of 
the regions in momentum space that dominate the low- 
energy properties of the material i.e. the Fermi surface. 
The expected large Fermi surface [as depicted in Fig. 
2(d)]has a volume greater than 1/2 (measured filling level 
<5 = 1.12 ± 0.05) and has a shape that is consistent with 
the LDA band calculations 11 1. Interestingly however, 
it shows two regions of suppressed spectral weight near 
(0.657T, 0.3-71") and (0.3-7T, 0.657r). Although there is obvi- 
ously some modulation of the intensity due to matrix el- 
ements, as some of the details of the intensity pattern are 
different between the two excitation energies, the gross 
features are the same. In addition, data from SSRL at 
16.5eV with the sample turned so that the polarization 
is at 45° to the Cu-0 bond reveals a similar pattern (not 
shown). The fact that the same systematics are seen in 
different configurations gives us confidence that we are 
measuring intrinsic properties, a fact which will be sub- 
stantiated by lineshape analysis below. 

A detailed look at Tllu — 16.5eV EDCs through the 
suppressed region of the Fermi surface [Fig. 1(d)] reveals 
that the peak initially approaches Ep and then mono- 
tonically loses weight despite the fact that its maximum 
never comes closer than ^100 meV to Ep. This results 
in features which are quite broad, even at kp. Such be- 
havior is similar to the high-energy pseudogap seen in 



the extreme underdoped p-type materials although in the 
present case it is observed near (0.657T, 0.3-7r) and not at 
(71-, 0) the maximum of the <i-wave functional form. An 
additional important difference is that a true clean gap at 
the very lowest energies cannot be defined, because the 
continually decreasing slope of the EDCs at kp gives a 
leading edge midpoint that is within a few meVs of Ep. 
We need to make the distinction between high-energy 
pseudogap behavior [the conversion of an expected sharp 
peak at Ep to a broad feature over a large energy range 
(~100 meV) and the suppression of low-energy spectral 
weight] and low-energy pseudogap behavior which is the 
clean leading edge gap seen in the normal state of under- 
doped and optimally doped p-types near (71", 0) that has 
been posited to be reflective of pairing fluctuations in the 
normal state. The low-energy pseudogap in the p-type 
materials appears to evolve somewhat continuously with 
underdoping into the high-energy pseudogap of the un- 
derdoped superconductor and AF insulator. Both gaps 
can be distinguished in moderately doped samples. The 
fact that they share the same d-wave functional form has 
lead to proposals that the pairing in the superconduc- 
tor is an inherited property of the insulating state [12| . 
These distinctions are important ones that we will return 
to later. 

In Fig. 3a and 3b we plot ?iw=55 and 16.5eV EDCs 
from around kp from (tt/2, it/ 2) to (it, O.Stt) as shown in 
the inset. Consistent with the intensity maps in Fig 2(a) 
and 2(b) the data show a well defined and intense peak 
at Ep in momentum space regions close to the T to (it, it) 
Fermi surface crossing and the (it, 0) to (71", it) Fermi sur- 
face crossing. In between these regions, corresponding to 
where the maps in Figs. 2(a) and 2(b) show suppressed 
intensity, the spectra are rather different [see curves 4 
of Fig. 3(a) and 7-9 of Fig. 3(c)] with no well defined 
peaks. These qualitatively different spectral lineshapes 
prove that the major effects seen in the intensity plots 
are not due to matrix elements which will not change the 
shape of the spectra drastically. These assertions receive 
further support from the measurement of EDC widths in 
Fig. 3(b). Here we plot the width in energy of the EDCs 
(defined as the FWHM/2 of the actual features and not 
the output of a fitting parameter) in Fig. 3(a). This 
quantity, which can be viewed as an average over the low 
energy scattering rate (ImE) is largest in the momentum 
region of interest along kp where the neai-Ep spectral 
weight is suppressed maximally as shown in Fig.'s 2(a) 
and 2(b). This analysis reveals those regions along the 
FS that show the anomalous low-energy spectral weight 
suppression are subject to the strongest scattering effects. 

We take the view that this anomalous electronic struc- 
ture in the suppressed region is reflective of the transfer of 
spectral weight from a near- i^p peak to a large broader 
incoherent maximum at higher energy. As one moves 
along the Fermi surface from (7r/2, it/ 2) to (it, 0.37r), this 
high-energy part at first gains intensity and shifts the 



global maxima to higher energy. In the intermediate re- 
gion along the kp contour, the near Ep weight is sup- 
pressed maximally with a plurality of spectral weight at 
higher energy. Moving towards (tt, 0.3tt) the high-energy 
part still forms a maximum at ~100 meV but the low- 
energy peak begins to recover at Ep and the two features 
can be seen simultaneously in a single spectrum. The 
low-energy peak continues to gain weight until the zone 
edge at ~(7r,0.37r), while the high-energy part loses in- 
tensity and disappears. The point of view that the two 
features are part of a single spectral function is supported 
by the fact that we observe only a single Fermi surface 
with the expected Luttinger volume. 

We note here that the regions of momentum space 
with the unusual low-energy behavior fall intriguingly 
close to the intersection of the underlying FS with the 
AFBZ boundary as shown by the dashed lines in Fig. 
2(a) and 2(b). This suppression of low-energy spectral 
weight and large scattering rate in certain regions on the 
FS bears resemblance to various theories that empha- 
size a coupling of charge carriers to a low-energy collec- 
tive mode. Central to these schemes is the supposition 
that there exists anomalous low energy scattering chan- 
nels with typically (tt, tt) momentum transfer. A simple 
phase space argument such as that shown in Fig. 2(d) 
shows that those charge carriers which lie at the intersec- 
tion of the FS with the AFBZ boundary will suffer the 
largest effect of anomalous (tt, tt) scatterings as these are 
the only FS locations that can have low-energy coupling 
with Q = (tt, tt). These heavily scattered regions of the 
FS have been referred to in the literature as "hot spots" . 
It has been suggested that the large back-scattering felt 
by charge carriers in the hot spots is the origin of the 
pseudogap in the underdoped hole- type materials p3| . 

A few candidates for such a scattering channel have 
been proposed. One realization would be a strong cou- 
pling of the charge carriers to (tt, tt) magnetic fluctua- 
tions. Indeed, strong AF fluctuations have been found 
at (tt, tt), as may have been expected in a material that 
is as close to the antiferromagnetic phase as this one j|. 
In this context the pseudogap is seen as a precursor to 
the formation of the antiferromagnetic phase. Other pos- 
sibilities include coupling to short range fluctuations of 
the CDW, DDW, or phononic type. A great deal of the- 
oretical effort has been devoted to schemes based on the 
above general considerations. [p~3| — 115|| - Of course we can- 
not rule out a combination of effects. Alternatively, large 
effective (tt, tt) scatterings can be caused by umklapp (U) 
processes which are only allowed for charge carriers on 
the AFBZ boundary. It has been suggested that these 
U processes create anomalously large effective backward 
scattering that leads to hot-spot phenomena and the cre- 
ation of a non- Fermi liquid without an additional broken 
symmetry. |lq ] 

In the underdoped p-type materials, due to band fill- 
ing considerations, the underlying Fermi surface's area is 



smaller and hence its intersection with the AFBZ bound- 
ary is closer to (tt, 0). This, along with the existence of 
the near Ep (tt, 0) extended saddle points, would cause 
these regions of the Brillouin zone to feel the largest ef- 
fect of (tt, tt) couplings if such processes exist. As alluded 
to earlier, this may lead to the formation of a pseudogap 
whose momentum dependence closely mimics and can 
be confused with the functional form of the <i-wave su- 
perconducting gap. In the present case of NCCO, the 
underlying Fermi surface is more distant from the (tt, 0) 
regions and the band there is flat 300 meV below Ep. 
In this scenario, the pseudogap regions would move away 
from (tt,0) on the FS and along kp towards (tt/2,tt/2), 
as observed. 

The class of theories that predict pseudogap forma- 
tion at the intersection of the AFBZ boundary with the 
FS contrast with those that explain high-energy pseudo- 
gap and hot spot effects to result from fluctuations of a 
very large energy superconducting pairing scale that has 
been inherited from the AF insulator. The latter class of 
theories predict that the momentum dependence of the 
pseudogap will follow the general trend of the d-wave su- 
perconducting gap i.e. attaining a maximum value at 
(tt, 0). The advantage of the present experiments on the 
electron-doped material is that the two different origins 
for pseudogaps give different functional forms, whereas 
they will give similar ones in the p-type materials as de- 
tailed above. In the very underdoped regime in the p- 
type materials, ARPES measures a high-energy d-wave 
like pseudogap on the order of 100 meV near the (tt, 0) 
region |17| ]. As we observe a high energy pseudogap in 
NCCO of a similar energy scale, but with a very dif- 
ferent momentum dependence, our results demonstrate 
the possibility that this very large d-wave-like pseudogap 
in the extreme underdoped p-type materials may not be 
related to preformed pairs, but instead may be due to 
pre-emergent magnetic order or other similar phenom- 
ena. Our observation still allows for the existence of a 
normal state low-energy pseudogap (defined by the lead- 
ing edge) that may be related to the existence of pairing 
fluctuations above T c . In this regard it appears that one 
must carefully discriminate between high-energy pseu- 
dogap behavior, which may be caused by the effects dis- 
cussed above and is tied to the intersection of the FS with 
the AFBZ boundary, and pairing fluctuations that may 
cause pseudogap behavior on a lower energy scale and 
should follow the d-wave functional form. If one wishes 
to interpret the physics of the p- and n-type cuprates in a 
comprehensive fashion then our results suggest that some 
of the confusing pseudogap phenomenology regarding en- 
ergy scales and doping dependencies may be reconciled 
by distinguishing effects of these kinds. 

We thank X.J. Zhou and Z. Hussain for beamline sup- 
port. The Stanford Synchrotron Radiation Laboratory 
and the Advanced Light Source are both operated by the 
DOE Office of Basic Energy Science, Division of Chem- 



ical Sciences and Material Sciences. Additional support 
comes from the Office of Naval Research: ONR Grants 
N00014-95-f-0760/N000f4-98-l-0f95. The Tokyo crys- 
tal growth work was supported in part by Grant-in- Aids 
for Scientific Research from the Ministry of Education, 
Science, Sports, and Culture, Japan, and the New En- 
ergy and Industrial Technology Development Organiza- 
tion of Japan (NEDO). The crystal growth at Stanford 
was supported by the U.S. Department of Energy un- 
der Contract Nos. DE-FG03-99ER45773-A001 and DE- 
AC03-76SF00515 and by a NSF CAREER Award. M.G. 
is also supported by the A. P. Sloan Foundation. 



[1] Y. Tokura, H. Takagi, and S. Uchida, Nature 337, 345- 

347 (1989) 
[2] F. Gollnik and M. Naito, Phys. Rev. B 58, 11734-11752 

(1998) 
[3] T.R. Thurston et. al, Phys. Rev. Lett. 65, 263-266 
(1990); K. Yamada et. al, Journal of Physics and Chem- 
istry of Solids 60, 1025-1029 (1999) 
[4] Z. Wang et. al. Phys. Rev. B, 43 3020-3025 (1991); S. 
Kubo, and M. Suzuki, Physica C, 185-189 1251-1252 
(1991); W. Jiang et. al, Phys. Rev. Lett. 73, 1291-1293 
(1994) 
[5] C.C. Tsuei and J.R. Kirtley, Phys. Rev. Lett. 85, 182-185 
(2000); R. Prozorov et. al, Phys. Rev. Lett. 85, 3700 
(2000); J. D. Kokales et. al, Phys. Rev. Lett., 85, 3696 
(2000) 
[6] N.P. Armitage et. al, Phys. Rev. Lett. 86, 1126 (2001) 
[7] Q. Huang et. al, Nature, 347 369-372 (1990) 
[8] Y. Onose et. al, Phys. Rev. Lett. 82 5120 (1999) 
[9] D. King et. al, Phys. Rev. Lett. 70, 3159-3162 (1993); 
R.O. Anderson et. al. Phys. Rev. Lett., 70, 3163-3166 
(1993) 
[10] A line of symmetrized data points was used immediately 
below and parallel to the zone diagonal as these points 
were missing from our initial data set. 
[11] S. Massida, N. Hamada, J. Yu, and A. Freeman, Physica 

C 157, 571-574 (1989) 
[12] G. Kotliar and J. Liu, Phys. Rev. B 38, 5142 (1988); H. 
Fukuyama, Prog. Theor. Phys. Suppl. 108, 287 (1992); 
P. A. Lee and X. G. Wen, Phys. Rev. Lett. 76, 503 (1996) 
[13] D. Pines, Z. Phys. B 103, 129 (1997); J. Schmalian, D. 
Pines and B. Stojkovic, Phys. Rev. Lett. 80, 3839 (1998); 
A.P. Kampf and J.R. Schrieffer, Phys. Rev. B 42, 7967 
(1990) 
[14] R. Hlubina and T.M. Rice, Phys. Rev. B 51, 9253 (1995) 
[15] S. Chakravarty et. al, Phys. Rev. B 63, 094503 (2001) 
[16] N. Furukawa, T.M. Rice, and M. Salmhofer, Phys. Rev. 
Lett. 81, 3195 (1998); G Honerkamp et. al, Phys. Rev. 
B 63, 035109 (2001) 
[17] D.S. Marshall et. al, Phys. Rev. Lett. 76, 4841 (1996) 





-0.8 -0.4 -0.6 -0.4 -0.2 

~\ 1 1 1 r, 





-0.8 -0.4 -0.8 -0.4 

Energy (eV) 

FIG. 1. Energy distribution curves from various directions 
in the Brillouin zone taken with hu) — 16.5eV at SSRL. (a) V - 
(n,n), (b)T- (tt,0), (c) (tt, 0)-(tt, tt), (d) (0.65tt, 0)-(0.65tt, tt). 
For a schematic refer to the black arrows in Fig. 2(a). 
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FIG. 2. (color), (a) and (b) Fermi surface of the partial 
Brillouin zone of NCCO taken with fkv — 16.5 and 55eV 
respectively. The plotted quantity is a 30 meV integration 
about Ef of each EDC plotted as a function of k. 16.5eV 
data were taken over a Brillouin zone octant and symmetrized 
across the F to (it, n) line, while the 55eV data was taken over 
a full quadrant [11]. The polarization direction is denoted by 
the double ended arrow. The dotted line is the antiferromag- 
netic Brillouin zone boundary, (c) LEED spectra of NCCO 
cleaved in situ at 10K. (d) Schematic showing only those re- 
gions of FS near the black circles can be coupled with a (n, n) 
scattering. 
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FIG. 3. (a) EDCs from along the Izf contour for huj=55 
eV. The graph plots (n/2, n/2) on the bottom and goes to 
(it, 0.371") at the top along the Uf contour given in the in- 
set, (b) AE (defined as FWHM/2) for EDCs in Fig. 3(a). 
(c) hui=16.5 EDCs along the same &f contour. The large 
momentum independent background [defined as the signal at 
(n, 7r)] has been subtracted out. 



